Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



An accretion model for the eclipsing cataclysmic variable 
PG 0859+415 



M. D. Still 



1,2 



1 Physics and Astronomy, University of St. Andrews, North Haugh, St. Andrews, Fife KY16 9SS (radsl@st-and.ac.uk) 

2 Astronomy Centre, School of Mathematical and Physical Sciences, University of Sussex, Falmer, Brighton BN1 9QH 



< 



> 
00 

o 

a: 
6 

C3 



Accepted 1996 May 20. Received 1996 April 25; in original form 1996 January 8 



ABSTRACT 

The emission lines found in the majority of cataclysmic variables are generally 
used as tracers of accretion flows which dominate the light at optical wavelengths. It 
has been suggested from previous observations that the eclipsing nova-like variable 
PG 0859+415 shows individualistic orbital characteristics which are inconsistent with 
canonical models of stellar accretion. We present spectrophotometry of this star which 
suggest that the standard picture is not in conflict with observation. We provide 
evidence that the shallow optical eclipses are of an extended bright spot rather than 
the accretion disc, and that the low-excitation lines are dominated by a transient 
absorption component, perhaps a result of the bright spot or accretion stream eclipsing 
the disc. We argue that this object may be another nova-like variable displaying long- 
term or cyclic variations in mass-transfer rate from the secondary star. 
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1 INTRODUCTION 

PG 0859+415 is a 14th magnitude, eclipsing nova-like vari- 
able with an orbital period of 3.7 h (Grauer et al. 1994). As 
such, it is a member of the cataclysmic variable family, con- 
sisting of a red dwarf (the secondary star) transferring mass 
onto a white dwarf companion (the primary star) via Roche 
lobe overflow. Since there is no observational evidence for a 
significant magnetic field about the compact object, accre- 
tion is assumed to occur via a disc, which must be stable to 
variations in mass transfer rate in order to explain the lack 
of observed dwarf nova outbursts. A recent observational 
review of the nova-likes is provided by Dhillon (1996). 

The majority of cataclysmic variables are emission line 
objects, which are used as tracers of the accretion flow that 
dominates the spectrum at optical wavelengths. In dwarf 
novae, emission features are superficially understood to be 
composites of lines from the accretion disc and the bright 
spot, where the gas stream from the secondary star impacts 
the outer rim of the disc (e.g. Marsh et al. 1990). Similar 
observations of the more energetic nova-likes, however, pro- 
vide a complex picture which is proving difficult to unravel. 
Optical lines appear to be composites of both emission and 
absorption line spectra from a host of possible sources, such 
as the disc (Rutten et al. 1994), a disc wind (Honeycutt, 
Schlegel & Kaitchuck 1986), the bright spot (Still, Dhillon 



& Jones 1995) , the gas stream overflowing the accretion disc 
(Hellier & Robinson 1994), and the irradiated inner face 
of the secondary star (Beuermann & Thomas 1990). Given 
such a range of possible contributions to optical spectra, it is 
unsurprising that nova-like emission shows such diversity in 
properties and that the nova-like sub-classification scheme 
is constantly being questioned. 

The blue excess object PG 0859+415 was confirmed as a 
cataclysmic variable during the optical campaign of Grauer 
et al. (1994; hereafter G94). Broad-band photometry showed 
rapid flickering and eclipses every 220 m, which was sub- 
sequently interpreted as the orbital period. Time-resolved 
spectroscopy reveals a high-excitation emission line spec- 
trum with the higher Balmer series containing shallow ab- 
sorption wings, presumably from the thick inner disc, and 
which display velocity modulation on the orbital period. The 
long-term constancy of light curves combined with the emis- 
sion profiles suggest a nova-like classification. 

Their data provided G94 with a number of provoca- 
tive ideas. Firstly, photometry shows that eclipses are pre- 
ceded by a broad hump, consistent with behaviour expected 
if the bright spot contributes to the optical flux. However, 
the depth of eclipse never exceeded the pre-hump photomet- 
ric level - suggesting the possibility that the eclipse is of the 
bright spot, rather than the disc. Secondly, study of Balmer 
line profiles reveal a prominent narrow emission source with 
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Table 1. Journal of observations. E is the cycle number plus 
binary phase with respect to the ephemeris derived in Sec. 3.1. 

Date Start End Start End No. of 

(1995 Feb) (UT) (E - 12 000) spectra 

16/17 23.66 3.90 326.43 327.58 46 
17/18 22.77 3.94 332.73 334.14 40 
18/19 22.84 3.98 339.29 340.69 35 

orbital-dependent velocities consistent with localized gas on 
the near side of the disc when observed at secondary star su- 
perior conjunction. Such a source is not predicted by models 
of disc accretion. 

In this paper we present spectrophotometry of 
PG 0859+415, consistent with the data obtained by G94, 
and suggest an acceptable emission model which is a com- 
posite of secondary star irradiation, gas stream absorption 
and emission from a wind and/or an accretion disc. We find 
further evidence that the eclipsed source is the bright spot 
rather than the disc, providing an opportunity to study vari- 
ations in mass transfer rate from the secondary star via 
changes in spot properties and disc size. 
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Figure 1. O — C diagram for the times of mid-eclipse of 
PG 0859+415. The dashed line defines the linear ephemeris re- 
fined in this paper, the square symbols denote the photometric 
timings of G94 and the triangles are the spectrophotometric tim- 
ings from the current data. 



2 OBSERVATIONS 

PG 0859+415 was observed on the nights beginning 1995 
Feb 16-18 with the 2.5 m Isaac Newton Telescope on La 
Palma. Exposure times began at 140 s on Feb 16 but were 
increased over the night to optimize signal and orbital res- 
olution. Dwell times of 400 s were decided upon and used 
for the rest of the run. Dead time for the dumping of data 
was generally 70 s between frames. The detector was a Tek- 
tronix CCD, windowed to 300 x 1124 pixels and mounted 
on the Intermediate Dispersion Spectrograph with a grating 
of 1200 lines mm" 1 . This configuration gave a wavelength 
range of approximately AA5811-6671 A at 1.7 A resolution. 
A journal of observations is given in Table 1. A nearby com- 
parison star was placed on the slit and several photometric 
spectra were obtained through a wide slit for the purpose of 
flux corrections. CuAr arc exposures were taken every 30- 
40 m in order to calibrate the wavelength scale and instru- 
mental flexure. The flux standard, Feige 66 (Massey et al. 
1988), was observed with a wide slit to correct for instrumen- 
tal response and convert the spectra to an absolute scale. 

Medium-scale sensitivity variations were removed with 
a balance frame prepared from tungsten lamp and sky 
flat-fields. Fourth-order polynomial fits to the sky were 
subtracted and raw spectra were then extracted using 
the optimal algorithm of Home (1986). Arc spectra for 
PG 0859+415, the comparison star and the flux standard 
were extracted from appropriate spatial locations on the arc 
frames. Fifth-order polynomial fits were made to the arc 
lines (rms scatter < 0.05A) and the wavelength scale for 
each spectrum was obtained by interpolating between two 
consecutive arc spectra. A comparison of the tabulated abso- 
lute flux values with a spline fit to the continuum of the flux 
standard provided a correction for large-scale instrumental 
response. Slit-loss corrections were achieved by dividing the 
PG 0859+415 spectra by a constant found from summing 
the flux under the continuum of the corresponding compar- 
ison star observation and multiplying by a constant derived 



Table 2. O-C timings for mid-eclipses of PG 0859+415, relative 
to the new ephemeris. 



Date 
(1995 Feb) 



HJD 



Cycle 
(E) 



O-C 

(<0 



17 2 449 765.5760(9) 12 327 -179(79) 

17 2 449 766.4943(3) 12 333 -50(24) 

18 2 449 766.6481(13) 12 334 28(110) 

19 2 449 767.5653(24) 12 340 63(213) 

similarly from the average of the comparison star wide-slit 
exposures. A wavelength-dependent correction proved to be 
impossible due to the low photon counts in individual com- 
parison star spectra. 



3 RESULTS 

3.1 Eclipse timings 

Four eclipses were covered by these observations, provid- 
ing an opportunity to test and refine the orbital ephemeris 
measured by G94. Eclipse profiles, which were produced by 
summing over line-free wavelengths of the continuum, were 
fitted with parabolae (details of light curve production may 
be found in Sec. 3.3). The minima of these, presented in Ta- 
ble 2, were combined with the photometric timings of G94 
and fit with a linear function, giving an apparently improved 
orbital ephemeris of: 

T = HJD 2447881.8584(7) + 0.15281249(2) E 

where To is the time of mid-eclipse, and E is the subsequent 
orbital cycle number. 

The residuals of this fit are presented in the O — C dia- 
gram of Fig. 1. It is unclear whether the presence of outliers 
is systematic or statistical, but their significance is discussed 
in Sec. 4. 
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Figure 2. The average spectrum of PG 0859+415. 



3.2 Average spectrum 

The average of all spectra of PG 0859+415 is presented in 
Fig. 2. Given that the radial velocity of either stellar compo- 
nent remains unknown, spectra have not been re-binned to 
remove orbital line smearing. We find a strong blue con- 
tinuum, superposed with weak emission lines of Ha and 
He I A5876 A. He I A5876 A displays an emission core with 
shallow absorption wings and the red wing is blended with 
interstellar absorption. The strong continuum and Hel pro- 
file suggest that a relatively large fraction of the disc is pro- 
jected on the sky, compared to the deeply eclipsing mem- 
bers of the class (e.g. UX UMa: Rutten et al. 1994). The 
Ha profile is single peaked, which is at odds with profiles of 
accretion disc emission, but this is not unexpected. 



3.3 Spectrophotometry 

A continuum light curve was produced by summing each 
exposure over the line-free wavelengths AA 5929-6257 A, 
AA 6358-6506 A and AA 6625-6646 A. This is presented in 
Fig. 3. Both the pre-eclipse hump and shallow eclipse are 
found in this light curve, although they appear less obvi- 
ous than in the magnitude plots presented by G94 because 
of the linear scale, the range chosen to present the data, 
and nightly variations. Consistent with G94's findings, the 
eclipse depth does not significantly exceed the pre-hump 
continuum level. 

Line light curves were produced by subtracting a spline 
fit to the line-free continuum from the data, and summing 
over the profiles. There is no significant evidence for eclipse 
above the noise limits of either feature. The Ha curve is 
similar to those of nova-likes whose lines are dominated by 
a component from the irradiated inner face of the secondary 
star (e.g. Dhillon, Jones & Marsh 1994). Given that stellar 
conjunction may not necessarily coincide with mid-eclipse 
(Sec 3.4), it is still true that maximum emission occurs at ap- 
proximately the superior conjunction of the secondary star 
when, the inner, irradiated face of the star is most visible. 
The same behaviour may be found in the He I A5876 A line, 
although the flux contribution is dominated by the absorp- 
tion wings, most probably originating in the inner accretion 
disc. 
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Figure 3. Orbital light curves for the continuum, Ha and 
HeiA5876A from PG 0859+415. Each curve is repeated over a 
second orbital cycle and orbital phase 0.0 corresponds to mid- 
cclipsc and is marked by an arrow; however, the dashed lines, 
which arc measured from Doppler images, are believed to be close 
to the true phase of stellar conjunction (Sec. 3.4). 



If the secondary star is the major component of line 
emission, it will be evident in the orbitally-resolved line pro- 
files. The data were averaged into phase bins, each covering 
1/10 of the orbital cycle. Multiple offsets were applied to 
each and the resulting stack is plotted in Fig. 4. There ap- 
pears to be at least two components to the Ha line, suggest- 
ing that a line model containing only secondary star emission 
is inadequate. The profile consists of two anti-phased peaks, 
one consistent with an origin from the secondary star, and 
the other consequently associated with the primary star. 
This second peak is the feature identified by G94 as emis- 
sion from a source at the front of the accretion disc, observed 
at superior conjunction of the secondary star. It should be 
noted however, that the profile can also be interpreted as 
one, or a composite of emission features, with a variable 
absorption component contaminating the core. Despite the 
small signal, the He I A4471 A feature displays the same be- 
haviour superposed over the broad absorption wings. 



3.4 Line profile fitting 

The data sampling and spectral resolution allow a more de- 
tailed analysis of the emission line profiles than the coarse 
phase bins of Fig. 4 allow. The data were continuum sub- 
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& Home 1988). This involves the mapping ol velocity-phase 
data (V,4>) to a velocity-velocity field (V x ,V y ) via: 
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Figure 4. Orbital phase-binned spectra of PG 0859+415, stacked 
with offsets which are multiples of 4 mjy. 



tracted in the usual way, re-binned onto a constant velocity 
scale of 39 kms -1 , and averaged into orbital phase bins each 
covering 2% of the orbit for Ha, and 5% for HeiA5876 A. 
The line profiles are presented as functions of orbital phase 
in Fig. 5. The linear grey-scale goes from white to black with 
increasing line intensity but has different zero points for each 
image, 0-3 mjy for Ha, and — 0.6-0. lmjy for HelA5876A. 
The data are repeated over a second orbital cycle and blank 
strips represent unfilled phase bins. 

We first consider the Ha feature, on which G94 base 
their analysis. The profile is dominated by two anti-phased 
emission components. The component with smaller velocity 
amplitude is approximately consistent with emission from 
the secondary star (Dhillon et al. 1994). Consequently, as- 
suming a stellar mass ratio (q = M s /M p ) less than unity, 
the anti-phased emission is consistent with a location on the 
back of the disc as viewed at the superior conjunction of the 
primary star, and as indicated by G94. Both emission com- 
ponents appear strongest at phase 0.5, as expected if the 
inferred spatial location of each source is correct. At orbital 
phases when the anti-phased emission does not dominate the 
profile, the line generally appears double-peaked, suggestive 
of a disc origin, or of absorption contamination from a lo- 
cal source within the accretion flow contaminating the core 
(Still et al. 1995). 

Further insight into the line sources is often provided 
by image reconstruction using Doppler tomography (Marsh 



f(v,4>)= / / i(v x ,v y ) X 

J — oc J — oo 

g(V--y + V x cos 4>-V y sin 4>)dV x dV y (1) 

where / is the line intensity at velocity V and orbital phase 
cf>, I is the emission distribution of the resulting map in ve- 
locity coordinates (V X ,V V ), 7 is the systemic velocity and 
g(V) is the local line profile. Examples of similar image re- 
constructions may be found in e.g. Marsh & Home (1990). 
Given that the profile of Ha possibly contains an absorption 
component, the preferred method is that of Fourier-filtered 
back- projection (Home 1992), thereby avoiding the positiv- 
ity constraint demanded on reconstructions by a maximum 
entropy approach. 

In order to reduce smoothing in the map, the data were 
convolved with a high-pass filter damped with a Gaussian 
cutoff to avoid high-frequency noise amplification. Back- 
projection results in the velocity map presented in Fig. 6. 
V x is the abscissa and V y the ordinate. Superposed over 
the image is an approximate binary configuration, where 
K s — 280 kms -1 and q = 0.5. The secondary Roche lobe 
is plotted, with the trajectory of the accretion stream, orig- 
inating at the inner Lagrangian point, and the Keplerian 
velocity about the white dwarf along the stream, originat- 
ing from close to the secondary star's centre of mass. This 
second trajectory represents the velocity of a Keplerian ac- 
cretion disc at the spatial location of the stream (Marsh 
et al. 1990). The three crosses represent the centres of mass 
for the secondary star, binary and primary star. An orbital 
solution for PG 0859+415 is not known with any confidence 
and the choice of binary parameters is discussed shortly. 

The features within this map are as expected from the 
trailed data. There are two major emission spots, one which 
we interpret as originating from the inner face of the the sec- 
ondary star and the other, if we consider the velocity field 
to be that of a Keplerian accretion disc about the white 
dwarf, from a local source on the back of the accretion disc. 
A weaker ring is approximately centred on the primary star, 
suggestive of accretion disc emission. Fits to the line pro- 
files, computed from the tomogram, are presented in Fig. 7. 
These are useful as a guide to the quality of the mapping 
procedure. Note, however, that pixel values in the map are 
not time- variable, and therefore the tomogram struggles to 
reconstruct the flux variations of individual line components 
within the data. 

Clearly, the secondary star emission in the Ha tomo- 
gram is not centred over the V x — axis. Provided we are 
willing to accept the secondary star as the origin for this 
feature, the most immediate explanation is that the eclipse 
ephemeris does not measure times of stellar conjunction. 
The consequent direction of coordinate rotation is consistent 
with the hypothesis of G94 that optical eclipses are of the 
bright spot only. A Gaussian fit to the secondary star peak 
indicates maximum emission at a velocity of 175 ± 4 kms -1 , 
and 0.064 ± 0.002 orbital cycles before mid-eclipse. Since the 
secondary star emission is most likely to originate on its in- 
ner face, the measured velocity may be taken as a lower 
limit to the stars radial velocity semi-amplitude. The mea- 
sured phase is an approximation of white dwarf superior 
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Figure 5. Trailed spectrograms of the Ho and HelA5876 A lines from PG 0859+415. The data has been repeated over a second orbital 
cycle and blank strips indicate unfilled phase bins. The two images are on separate linear intensity scales. 
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Figure 6. Dopplcr tomograms of the Ho and HclA5876A feature from PG 0859+415. The intensity scales are linear but different for 
the two maps, where white to black illustrates increasing flux. 



conjunction, which is illustrated by both the orientation of 
the binary on the tomograms and the dashed line found in 
Fig. 3. However, if irradiation of the secondary star is not 
symmetric about the Li-point, perhaps a consequence of 
shadowing by the bright spot or accretion stream (South- 
well et al. 1995), then the conjunction phase will be inac- 
curate by no more than a few percent of the orbital period. 
The Ha data are not suitable for providing a measurement 
of the primary stars radial velocity semi-amplitude since it 
remains unclear whether any of the emission is distributed 
symmetrically about the white dwarf (Shatter 1985; Still 
1996). Therefore the binary parameters adopted for Fig. 6 
are a guess, based on a mass ratio of q = 0.5 and assuming 
that only the inner face of the secondary star is irradiated. 

The emission spot at the back of the disc has been in- 
ferred from two independent data sets, yet there appears no 
plausible reason for its presence. A solution is provided from 
an inspection of the He I A5876 A profiles, which suggest that 
the spot is an artifact of core absorption. Fig. 5 illustrates 
that the two emission features found in the Ha line are also 
observed in the Hel profile. Also present, however, is an ab- 
sorption component which dominates the line core and has 
variable velocity and intensity over the orbital cycle, attain- 
ing maximum strength between orbital phases 0.6 and 1.1 
and maximum velocity at phase 0.9. This provides a con- 
sistent, yet alternative, model for the line distribution: The 
first component is a broad emission feature which is visi- 
ble in Fig. 5 between —500 kms -1 and +500 kms -1 ; this 
is most likely from the accretion disc (Rutten et al. 1994), 



or an accretion disc wind (Dhillon, Marsh & Jones 1991). 
The second is an absorption line overlaying the first compo- 
nent, similar to the feature observed in the deeply eclipsing 
nova- like RW Tri (Still et al. 1995), and the third is nar- 
row emission from the secondary star (Dhillon et al. 1994). 
The last component appears to be faint but still visible at 
phase 0.9 where it is observed to cross the absorption fea- 
ture. This is plausible given the inclination inferred by the 
shallow eclipse, and the superposition of this narrow com- 
ponent over the absorption provides confidence that it is an 
individual emission component and not the residue of core 
absorption. 

With hindsight, this line model is equally plausible for 
the Ha feature. The inferred spot at the back of the disc 
is a consequence of the same absorption component, found 
weakly in Ha, which depletes the broad emission component, 
leaving a residue which maps to the extended spot in Fig. 6. 

Before back-projecting the suitably-filtered 
He I A5876 A data, the interstellar features contaminating 
the red wing of the line were interpolated over. The result- 
ing map, presented in Fig. 6, shows the features interpreted 
in the data as the secondary star component, the strong 
absorption core and the broad emission residual. Since the 
individual pixels within the map are time-independent, the 
narrow secondary star emission and the absorption core have 
not been reconstructed accurately. Hence, as is common, for 
nova-like emission, the computed profiles of Fig. 7 do not 
provide a good representation of flux variations in the data. 

If incorporated into the model of Keplerian disc and 
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ballistic accretion stream, both the map and computed pro- 
files indicate that the absorption occurs on the trailing side 
of the disc, although it does not appear to coincide with ei- 
ther the gas stream trajectory or the Keplerian trajectory 
off the stream. Although there is some freedom to arrange 
the conjunction ephemeris and binary parameters accord- 
ing to the irradiation distribution on the inner face of the 
secondary star, there are no plausible configurations to al- 
low the stream trajectories to coincide with the absorption 
velocities. However, absorption strength appears correlated 
with the continuum hump displayed by the optical contin- 
uum (G94 and Fig. 3), suggesting that the feature could be 
related to the bright spot or accretion stream. A possible 
interpretation of its velocity is that the disc is thicker down- 
stream from the bright spot than upstream, and the disc 
rim at these azimuths is able to eclipse the inner disc - con- 
sistent with the model suggesting that it is the bright spot 
which is being eclipsed by the secondary star. The absence 
of absorption at the likely bright spot velocity is perhaps 
a consequence of core filling by emission from the spot or 
stream. The data is not of the desired quality to favour or 
refute this idea; more precise light curves of the lines features 
during eclipse are required for this purpose. 

Certainly, the line distribution does not provide a 
unique emissivity model, however, the strength of the above 
interpretation is that it considers only accretion elements 
that are well-founded in both theory and observation. 



4 DISCUSSION 

The present data can be described by a model where the 
bright spot is the major contributor driving orbital phe- 
nomena in the nova-like variable PC 0859+415. The bright 
spot appears to be eclipsed by the secondary star, and low- 
cxcitation lines are dominated by a variable component, 
perhaps from the spot region, which occults the inner disc 
and provides a confusing set of emission profiles to decon- 
volve. Fig. 1 illustrates that there are a number of outlying 
eclipse times from the best linear orbital ephemeris. It is 
not implausible that these outliers are the result of mislead- 
ing statistics. However, it may be that this distribution is 
not best fit by a linear ephemeris and is a result of variable 
mass transfer from the secondary star (Robinson, Shetrone 
& Africano 1991; Applegate 1992). If the emission model 
suggested by this data proves to be accurate, long-term pho- 
tometric and spectroscopic monitoring of bright spot inten- 
sities and relative disc sizes, provided by eclipse times, may 
provide an opportunity to study the accretion flows response 
to variations in mass transfer rate from the secondary star. 
Further eclipse times are required to improve the statistics 
of the O — C diagram before the significance of the linear 
ephemeris may be suitably tested. 



5 CONCLUSIONS 

The UX UMa nova-like classification made by G94 is 
consistent with the observation made in this paper that 
PG 0859+415 is a lower inclination analogue of UX UMa 
(Kaitchuck, Honeycutt & Schlegel 1983) and RW Tri (Still 
et al. 1995). Lines are interpreted as composites of the bright 



spot and inner accretion disc in absorption, and the irradi- 
ated inner face of the secondary star and outer disc, or a 
disc wind, in emission. If unrecognised, the bright spot com- 
ponent, provides a confusing emission distribution which is 
at odds with models of disc accretion. Provided the identifi- 
cation of secondary star emission is correct, a lower limit to 
the radial velocity of the donor star is 280 kms -1 , and the 
eclipse found in the optical continuum must be dominated by 
the bright spot rather than the disc. Further eclipse times 
are required to determine the significance of outliers from 
the best linear orbital ephemeris, which may be indicative 
of variable mass transfer from the secondary star. 
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